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Aspects of Density Functional Resonance Theory (DFRT) [Phys. Rev. Lett. 107, 
Ph 163002 (2011)], a recently developed complex-scaled version of ground-state Density 

Functional Theory (DFT), are studied in detail. The asymptotic behavior of the 
complex density function is related to the complex resonance energy and system's 
Sp-i threshold energy, and the function's local oscillatory behavior is connected with pref- 

erential directions of electron decay. Practical considerations for implementation of 
the theory are addressed including sensitivity to the complex-scaling parameter, 9. 
In Kohn-Sham DFRT, it is shown that almost all ^-dependence in the calculated 
> energies and lifetimes can be extinguished via use of a proper basis set or fine grid, 

i The highest occupied Kohn-Sham orbital energy and lifetime are related to a physical 

affinity and width, and the threshold energy of the Kohn-Sham system is shown to 
be equal to the threshold energy of the interacting system shifted by a well-defined 
y—( functional. Finally, various complex-scaling conditions are derived which relate the 
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functionals of ground-state DFT to those of DFRT via proper scaling factors and a 
^ non-Hermitian coupling constant system. 
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I. INTRODUCTION 



Density Functional Theory (DFT) 1-3 has provided one of the most accurate and reliable 
methods to calculate the electronic properties of molecules, clusters, and materials from 
first principles. It is the most widely used formalism in modern computational quantum 
chemistry 4 . In addition, DFT's time-dependent extension (TDDFT) 5 can now be applied 
to a wealth of excited-state and time-dependent properties in both linear and non-linear 
regimes 6 . An advantage of DFT over wavefunction-based methods is that the one-body 
electron density is the primary variable, not the intracable many-body wavefunction. 

Recently, a complex-scaled version of DFT, Density Functional Resonance Theory 
(DFRT), was introduced 7 that allows one to calculate the in-principle exact resonance 
energy and lifetime of the Lowest Energy Resonance (LER) of an unbound system. This 
theory makes use of a complex "Kohn-Sham" system that facilitates self-consistent calcu- 
lations on many-electron systems with a complex-valued density function as the primary 
variable. As shown in Refs. 7,8 and in this work, the use of a localized complex density 
is advantageous, especially in the case of resonances where the wavefunction is not only 
an awkward many-body object but also divergent, or non-normalizeable, because of the 
transient nature of the system. As we will argue here, DFRT is a promising approach to 
calculate negative electron affinities (NEA) and resonance energies and lifetimes. 

Shape and Feshbach resonances in low-energy electron scattering processes 9-11 are of grow- 
ing interest in biological systems 12-14 , atmospheric sciences, lasers, and astrophysics 15-18 . 
For example, consider electron impact on DNA. Experimentalists have made great strides 
in their understanding of damage to DNA caused by secondary low-energy electrons during 
cancer radiation. It turns out that this secondary process is much more important than 
previously thought 19 . Next to these experimental efforts, theoretical calculations are either 
very limited or extremely non-trivial due to the size of the relevant species, the many-body 
nature of the problem, and the bound-free correlations. Even in small systems such as the 
molecular hydrogen anion or doubly negative ions, stability and structure are still a very 
active subject of research 10,20 . 

In addition to relevant metastable properties that are within reach of the theory, DFRT 
reduces to standard ground-state DFT with the removal of the complex-scaling transfor- 
mation and thus allows one to examine current topics of research in approximate DFT 
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from a different viewpoint. One example of such an application is recent work investigating 
derivative discontinuities of the energy at the maximum number of bound electrons 21 . 

In this paper, several aspects of the DFRT formalism, introduced in Ref. 7, are studied 
in detail. The outline of the paper is as follows. First, in section II we review the analog of 
Hohenberg and Kohn's theorems for the LER and the main ideas of DFRT. Aspects of con- 
vergence in DFRT calculations are addressed in section III. Section IV provides a discussion 
of the physical significance and use of the complex density function, the primary variable 
of DFRT. The Kohn-Sham orbital energies and lifetimes are related to physical properties 
of the interacting system in section V, and possible approximations for the complex uni- 
versal energy functional are explored in section VI. Lastly, we conclude and mention some 
promising applications and extensions of DFRT. 



II. DENSITY FUNCTIONAL RESONANCE THEORY 

To begin, we give a brief sketch of the analogs of the Hohenberg-Kohn (HK) theorems 22 
for a LER, and we describe the main ideas of Kohn-Sham DFRT. 

A quantum resonance is associated with a complex number, E n = S n — (T n /2)i, that is 
a pole of the scattering matrix or a peak in the continuum density of states. The index n 
labels a specific resonance (atomic units are used throughout). £ n is the resonance energy 
or position, and T n is the resonance width. The lifetimes of the resonances, C n , are given 
by the inverse of the widths. Consider the LER of an unbound system and assume that 
this lowest energy metastable state is also the longest-lived resonance (this condition is not 
proven, but is known to be the typical case 23 ' 24 ). 

The complex "density" associated with the LER is 

n,(r) = <*£|n(r)|*?> (1) 

where n(r) is the density operator, and (^| and \^$) are the left and right eigenvectors of 
the complex scaled Hamiltonian, Hg, corresponding to the LER. The angle 9 is the complex 
scaling parameter in the transformation of the coordinates from r to re te (see Ref. 25 for a 
review of such transformations). 

In general, there is no variational principle for the complex energy associated with the 
LER, but only a stability principle. In Ref. 22, a complex density variational principle was 
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presented for the LER pertaining to trial functions of a certain kind. Specifically those trial 
functions, <&triah which can be expanded in a set of resonance wavefunctions and for which 
\(^e\^f r i a i)\ > 1/2- In addition, Ernzerhof and co-workers have shown that for potentials 
and densities with a small imaginary part (compared to their real part), the standard density 
variational principle is maintained 26 . 

A one-to-one correspondence between complex densities and complex-scaled potentials 
can be shown via contradiction, but the Levy-Lieb constrained search approach is taken 
here. In this approach, the lowest-energy resonance energy and lifetime of a system defined 
by the Hamiltonian, H — T + V ee + J drn(r)v(r), can be written as: 

E hEK = min 

x(^\f e + V e e e + J dvn{v)v{ve te )\^) (2) 

where the energy and width are minimized over all resonance wavefunctions. Next, split 
the minimization into two steps. First, search over all resonance wavefunctions that give a 
certain complex density, and then search over all complex densities. 
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Here, F d [no] is a complex- valued "universal functional" defined by, 
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F°[n e ] = min I | (^|T» + KelO (4) 

Ve^n e I _ 2Im 

Carrying out the minimization with respect to variations in the density given the constraint 
that the density integrates to the number of electrons (/i - Lagrange multiplier) gives 
5 
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and one obtains the result, 
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Therefore, v (re ie ) maps directly to the specific n e (r) that satisfies the above condition. 
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In the original extension of the Hohenberg-Kohn Theorems to LER's 22 , both a proof by 
contradiction and the Levy-Lieb type constained search were employed. In the proof by 
contradiction, it is assumed that the energy and width are minimized over all trial complex 
densities that are v e -representable. This means that the trial densities can be associated with 
a Hamiltonian containing the smooth complex-scaled external potential v(re td ). However, if 
an arbitrary trial complex density is chosen, it is not immediately evident that this condition 
is fulfilled. The Levy-Lieb constrained search of Eq. 3 lifts this restriction and replaces it 
with a weaker restriction on the trial complex densities. The search in Eq. 3 is performed 
over all iV-representable densities. This is very similar to the iV-representability condition 
normally applied in ground-state DFT (see Ref. 27 for a review of f-representability and 
7V-representability in ground-state DFT), but here we search over all localized complex- 
valued functions that integrate to the number of particles, rather than searching over purely 
real functions. This is a larger space of densities than what is taken in ground-state DFT 
and hence the iV-representability condition of DFRT is a weaker restriction than that of 
iV-representability in DFT. 

Having motivated this correspondence between complex densities and complex-scaled 
potentials, we proceed with defining the theory. It is required that n$(r), corresponding to 
the LER, be normalized to the number of electrons, as real densities are: 



The energy and lifetime of the resonance can be extracted from rig with a properly scaled 
energy functional. For N electrons we write this functional as 



where v is the external potential, T s e the complex-scaled non-interacting Kinetic energy 
functional, is the classical Hartree energy functional, and E xc is the exchange-correlation 
energy functional. We require that T®[ng] = e~ 2ld T s [ng\ and ^[^e] = e E H [ng], where 
T s [ne] and E H [ng] are the standard non-interacting kinetic energy and Hartree functionals 
evaluated at the complex densities. Without an explicit expression for E^ c [ng], however, the 
total energy cannot be calculated via Eq. 8. The exact scaling of functionals and possible 
candidates for approximate functionals will be discussed in Section VI. 




(7) 




+E e H [n e }+E e xc [ne] 



(8) 
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Now map the system of interacting electrons whose LER density is ne(r) to one of 
N particles moving independently in a complex "Kohn-Sham" potential v e s (r) defined 
such that its N occupied orbitals {4> 9 (r)} yield the interacting LER-density via ng(r) = 
YliLii&t' L |^( r )l0f R )- ^ n Moiseyev's Hermitian representation of complex-scaling 28 , the 
complex Kohn-Sham equations are: 



hi - Si -h 2 - 2T7 1 \ I Re(0f ' e ) 
h 2 + 2rr l h x -e i ) \ lm(0f 



, (9) 



where h x = -\ cos(2#)V 2 + Re(^(r)), and h 2 = \ sin(2#)V 2 + Im(t;f(r)). The set of and 
{-Tj} are the orbital resonance energies and lifetimes of the Kohn-Sham particles, and the set 
of {4>f' 6 } and {(f)f' 9 } are the right and left eigenfunctions of the non-Hermitian Kohn-Sham 
Hamiltonian. 

To show that the solution of Eq. 8 can be written in an Euler-Lagrange form, expand the 
right and left Kohn-Sham orbitals in an orthonormal basis, 

N N 

i i 

The complex density can then be written as, 

N 

ne(r) = J2cfc*xKr) (H) 

i 

and the functional derivative of E e [ng] = £[ng] — ^C~ l [ne] with respect to n$(v) has the form 

5E 9 [n e } 1 dE e 1 8E e 

5n e (r) ~ 2cfx?( r ) dc i 2cfx?( r ) dc i 
1 dE e 

Moiseyev et al. have shown that the resonance wavefunction is given by setting 8E e /dcf = 
dE 9 /dcf = for all i 25 . Therefore, the LER solution of Eq. 8, with the constraint that the 
complex density integrates to the number of particles, is given by, 

SE e [n e ] 



J drn e (r) = . (13) 



5n e 

Performing the variarion in Eq. 8 and comparing with Eq. 9 leads to an expression for the 
Kohn-Sham potential that is again analogous to that of standard KS-DFT: 

v e s (v) = v(ve id ) + e-»v*[no](T) + vi c [n e ](v) , (14) 



where v^ c [ng](r) = 5E^ c [n g ]/Sng(r)\i I ER. Now with this auxiliary Kohn-Sham system estab- 
lished, self-consistent calculations for the LER's energy and lifetime are possible with an 
intial guess for the complex density. 

III. ASPECTS OF CONVERGENCE IN DFRT CALCULATIONS 

The main structure of the DFRT formalism mimics that of ground-state DFT. One 
chooses an initial guess for the complex density, calculates the KS potential from Eq. 14, 
solves the KS equations of Eq. 9, constructs a new density by summing the KS orbitals, and 
repeats until self-consistency is reached. Despite the similarities, Kohn-Sham DFRT is a 
complex-scaled Non-Hermitian formalism, and a number of unique practical aspects of the 
theory should be discussed. 

The resonance energy and lifetime is independent of 9 in complex-scaling theory 29 . How- 
ever when doing a numerical calculation, the use of finite basis sets or finite grids necessitates 
finding an optimum 9. This optimum condition is acheived by examining ^-trajectories 29 . 
When the real and imaginary parts of the resonance energy become stationary around the 
optimum value of 9 one can see a kink or loop in the trajectory. With many more grid 
points or a larger basis this procedure becomes less important as all the ^-trajectory points 
collapse to the region of the resonance. 

^-independence of the energy is preserved by the SCF procedure of DFRT. As the 
grid or basis size increases the dependence on 9 becomes negligible. To illustrate this 9- 
independence, calculations were performed on a model potential similar to the one discussed 
in Ref. 7 filled with two solf-coulomb interacting electrons. The potential has the form, 



where a, a, 6, c, and d are constants. Results for the energy and lifetime of a two-electron 
resonance were obtained using the Fourier Grid Hamiltonian (FGH) method and with both 
a particle-in-a-box and harmonic oscillator basis. It can be seen that the method of using 
rig in a scaled functional is independent of 9, as long as the basis is large enough. In other 
words, for a given choice of 9 the density obtained using a large grid and a large basis will be 
the same, and those densities will yield the same result in the scaled functional. The error 




(15) 
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FIG. 1. The difference between the real and imaginary parts of ng from the basis function and 
the fourier grid method. Given a choice of 6 these converge to the same density, (a = 4, b = 0.5, 
c = 4, and d = 2) 

in the density, plotted in Fig. 1 versus the number of basis functions, was calculated as, 

error Ro = J dx (Re(n P ) lb (x)) — Re(riQ gh (x))) 2 

error Im = J dx {lm{n p e lb {x)) — \m.{nl gh {x))) 2 (16) 

where rJ e 9h is an accurate density from the Fourier grid Hamiltonian method and is the 
density obtained using the particle-in-a-box basis (both obtained with the same 9). Figure 1 
shows the convergence of the density as the number of basis functions increases. Note that 
if one chose a smaller basis or grid, an optimum 9 could still be found for the calculation. 
This optimum 9 gives the correct resonance energy and lifetime 29-32 . This principle of 
^-independence is important, because within a Kohn-Sham DFRT calculation one has to 
make a choice for the scaling parameter. Yet, the advantage of the scheme is that one only 
has to solve one-body, non-interacting equations. For such equations, one should be able 
to efficiently use a large enough basis set or a fine enough grid to extinguish most of the 
numerical ^-dependence. Thus, a well-known drawback of the complex-scaling technique 29-31 
is outdone by the benefit of never having to deal with iV-particle wavefunctions, but just 
one-body (complex) densities. 

Although the energy in DFRT is theoretically independent of the scaling parameter, all 
of the potentials and the density itself depend on the choice of 9. Therefore, one should be 
careful to compare only potentials and densities calculated with the same 9. 



In addition to the issue of ^-dependence, convergence of the solutions should be discussed. 
As can be seen in Fig. 1 and in Ref. 8, often the imaginary part of the energy has a slower 
convergence than the real part. This occurs when the width (Im(E e )) of the resonance is 
much smaller than the position (Re(E )) which is typically true for long-lived resonances. 
One must evaluate the functional using a very accurate complex density in order to precisely 
calculate a narrow width. However, the self-consistent DFRT method has already been 
tested for small model systems of interacting electrons where this accuracy is possible 7 ' 8 . 



IV. THE COMPLEX DENSITY FUNCTION 

A natural question is: what kind of physical significance can be associated with the 
real and imaginary parts of the complex density function? This question has already been 
partially explored by Moiseyev and Barkay. They showed that the phase of the square 
root of the complex density in systems like quantum dots or tunneling diodes is related 
to the phase of a measureable quantity, the complex tunneling probability amplitude 33 . 
This quantity provides the probability for obtaining a specific scattering result such as 
transmission through a diode. Using this property of the complex density, Moiseyev and 
Barkay derived a formula for the tunneling probability amplitude which correctly predicted 
the energetic position of resonances. In addition to the study by Moiseyev and Barkay, 
Buchleitner et al have examined how a properly normalized magnitude of the complex 
density can be interpreted as the real electron density of the metastable system 25 ' 34 . The 
fact that the real part of the density might resemble a bound density for very long-lived 
states has been emphasized in previous work 8 . 

In addition to these interpretations, the asymptotics of the LER's complex density contain 
relevant information about the real system. As r — > oo the density, ne, behaves like 25 

Mr) ~ Ae l2 ^ eWr (17) 

Where AE = E 6 — E th , E 9 is the complex energy of the LER, and E th is the threshold 
energy for the relevant decay channel. The assymptotic expression for ne can be split up 
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into real and imaginary parts. 

n e (r) ~ Ae -^\me*°^E)r 



Re 
Im 



, cos2 v / 2Re(e if? v / AE)r , 

(18) 

sin 2 v / 2Re(e* e v / A£)r 

One can see that the decay of the real part of rig is coupled to both the real and imaginary 
parts of AE as is the imaginary part of rig. If the density is written in the following way 

n e (r) = \n e (r)\e^ (19) 

where at large r, 

\ng(r)\^Ae- 2 ^ lm ^ e ^> (20) 

0(r) ~ 2v / 2Re(e ie v / AE)r (21) 

then the decay of the magnitude of rig is governed by only the imaginary part of e t0 AE and 
the decay of the phase of rig is governed by only the real part of e t6 AE. 

Along with the magnitude and phase of rig, the structure of the complex function itself can 
provide physical insight. To illustrate this point consider soft-coulomb interacting particles 
in a ID potential, v(x) + jQ(—x), where 7 is a constant and v(x) is given in Eq. 15. This 
potential only supports resonances and has a bias for the decay from one side to another 
depending on the value of 7 (see Fig. 2). The complex density function was calculated, and 
an example density is included in Figure 3. Strong oscillations in the tails of the density 
indicate the geometrical bias of decay. In this simple example, an electron prefers to exit 
the potential to the right. Therefore, this system will tend to interact, or react, more with 
another system that is brought in from the right rather than from the left. This structure 
motivates the use of the complex density function in studying the reactivity of negative ions. 
In three dimensions a gradient of the complex density could be employed to see local areas 
where there is a geometrical bias of decay. It is out of these critical areas that the metastable 
system would tend to donate an electron in a chemical reaction. Metastable systems are 
often very reactive due to their diffuse electron clouds, and a systematic theoretical method 
for exploring reactivity would be of general interest. This is related to work from Moiseyev 
showing how one can extract information about the partial widths from the decay of the 
complex wavefunctions in assymetric potentials 35,36 . The advantage here is the use of the 
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FIG. 3. Real and imaginary one electron complex densities. One with an assymetric bias and one 
without (a = 7, b = 0.05, c = 4, d = 1.5). 

one-body complex density and not the many-body resonance wavefunctions. It is possible 
to develop a useful reactivity theory based on DFRT calculations for the complex density 
and its various derivatives. 



V. KOHN-SHAM ORBITAL ENERGIES AND LIFETIMES 

As in DFT, the auxilary KS orbitals or energies in DFRT are assumed to not have a direct 
physical interpretation. However, in recent work, the authors have shown that the highest 
occupied orbital energy and lifetime of DFRT can be related to physical quantities 21 (similar 
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to Koopmans' theorem of DFT). Specifically if the system of interest is a (N + l)-electron 
metastable system with a single dominant decay channel, 

(4-^) = (-^-^) , (22) 

where e e H = e H — 2ir H 1 is the complex Highest Occupied (HOMO) orbital energy of DFRT, 
En is the HOMO resonance energy or position, r H is the HOMO resonance lifetime, A is the 
negative electron affinity (NEA) of the iV-electron system, T is the width of the metastable 
system, and e th is the KS "threshold" energy defined by, 

e th = E th _ ^ ( 23 ) 

where £[ng] = J^Ii e f + ^Hxc[ n e] ~~ f drv^ xc (r)ng(r). In other words, the threshold energy 
in the KS system is the threshold energy in the real system shifted by — £. Individually the 
KS HOMO resonance energy is Re(4) = e H = Re(-A) + Re(e th ) and the KS HOMO width 
(inverse lifetime) is Im(e^) = 2-r" 1 = — | + lm(e th ). Both the energy and lifetime are real 
physical quantities shifted by either the real or imaginary part of e th , a factor resulting from 
Hartree, exchange, and correlation contributions. 

For a bound system, the standard Koopmans' theorem of DFT is recovered and the 
HOMO energy of DFRT is equal to the negative of the ionization potential of the system. 

We emphasize that DFRT is applicable to both bound ground-states and LER's. There- 
fore, a system of interest might have some KS orbital energies that are bound and some 
that are resonances. For example, consider a metastable negative ion formed by adding 
an electron to a neutral bound system having a NEA. A DFRT calculation on the bound 
parent system would render all orbital energies to have zero imaginary part (i.e. all the 
Kohn-Sham particles are bound). A DFRT calculation on the metastable anion would yield 
all but one KS particle bound. The remaining particle (the HOMO) would be a resonance 
with its resonance position and width given by the expressions above. 

In addition to relating the KS HOMO energy to physical quantities, a useful feature of 
DFT that is preserved in DFRT is the ability to construct the total energy as a corrected 
sum of orbital energies, 

N 

E e [n e ] = Y t (^-^r 1 )+E e HX [n 9 ] 
i=i 

-J drvi x (r)n e (v) (24) 
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This is a useful feature for practical implementations. 

A peculiarity of the Kohn-Sham lifetimes can be seen by considering a system of non- 
interacting particles in a potential that only supports resonances. In this case, the Kohn- 
Sham potential of DFRT is just the complex-scaled external potential, and Eq. 24 simply 
reduces to the sum over orbital energies. The lifetime of the full system is C^ 1 = J^=i r i _1 - 
For two non-interacting particles this implies that C = r H /2 and the lifetime of the system 
of two particles is less than that of the system with one particle even though the particles 
are non-interacting. 

VI. SCALING OF EXACT AND APPROXIMATE FUNCTIONALS 

Any implementation of DFRT must approximate the universal functional of Eq. 4. Re- 
lated work by Ernzerhof 26 and physical intuition suggest that bound ground-state functionals 
are applicable. If the energy functional is partitioned as in Eq. 8, the complex-coordinate 
scaling of the kinetic and Hartree functionals follows immediately from the coordinate scal- 
ing of the kinetic and interaction operators: T^[n e ] = e~ 2%e T s [uq] and E^[n e ] = e~ t8 E H [ng], 
where T s [ng] and £^[^0] are the standard non-interacting kinetic energy and Hartree DFT 
functionals evaluated at the complex densities. The scaling of the exchange and correlation 
functional is not immediately apparent. 

To study the complex-scaling of functionals, one can take advantage of the fact that 
DFRT can be applied to bound ground-states, not just LER's. For a system with a bound 
ground state, various coordinate scaling relations (r — >■ ar for a > 0) are known for the 
functionals of ground-state DFT 37 . Our strategy to derive exact relations between DFT and 
DFRT functionals, is to consider the complex-coordinate-scaling transformation, r — > re ie , 
keeping in mind that DFT and DFRT give the exact same energy for a system with a bound 
ground-state. For such a system, this correspondence between the energies of DFT and 
DFRT implies: 



where F[n] and Fo[ne] are the universal functionals of DFT and DFRT respectively, and 
n(r) and n e (r) are the real- valued density of DFT and the complex- valued density of DFRT 
corresponding to the ground-state. Once the complex-scaling relations of the functionals are 
in place, the assumption is made that the scaling does not change for a LER, even though 




(25) 
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Eq. 25 is no longer valid. This assumption is reasonable because the functional (Eq. 4) of 
DFRT is applicable to both bound and LER states. 

Now, consider the change in the ground-state functionals with a complex-scaling of the 
wavefunction and density coordinates. As in standard DFT, it is required that the scaled 
density integrates to the number of electrons. 

J dr n e (r) = J dr n(r) = N (26) 

In three dimensions this means: 

n(r) -> n e {r) = e 3ie n(re ie ) (27) 

so that the transformation (coordinate scaling) preserves the number of electrons. For 
the non-interacting kinetic energy, Tg[ng], the ground-state wavefunction that yields Uq 
according to Eq. 26 is given by 

, (n, ...,i*) = e 3im ^(v ie i9 , r N e w ) (28) 

Then, 

(V e \f \*o) = e 2ie (^|f |*) (29) 

where the "c-product," indicated by (^ e \^/ e ), is used to emphasize that the imaginary part of 
tyg that becomes complex purely from the transformation is not conjugated 25 . Eq. 29 shows 
that the functional has a coefficient of e 2ie from the scaling of the wavefunction. Yet, the 
complex-coordinate scaling leaves bound state energies untouched. Therefore, the scaling 
factor in the complex-scaled functional of DFRT, which is also applicable to bound states, 
must compensate for this factor of e 2td , 

T e s [n e ] = e^ e T s [ng} (30) 

Next, consider the Hartree functional. Like the kinetic functional, the relation of the complex 
Hartree functional of DFRT, E^[n e ], to the standard ground-state Hartree functional, £^[71], 
follows from the coordinate scaling of the interaction operator. To emphasize this point, the 
change in E n [n] caused by a scaling of the density coordinates is considered. By directly 
substituting n e (r), defined in Eq. 27, into E^bi], 

E s [n ]=e»E H [n] (31) 
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and the functional must compensate for this factor, 

E t\ n e] = e~ ie E H [n g ] (32) 

The exchange functional, E x [n], is defined by, 

E x [n] = (*'\V ee \V a )-E s [n] (33) 

where ^ s is the Kohn-Sham wavefunction that yields n(r) and minimizes the expectation 
value of T. Using Eq. 31 together with Eq. 28, 

E x [n e ] = e i9 E x [n] (34) 

and, 

Ei[n e ] = e- ld E x [n e ] (35) 
Lastly, consider the correlation functional, E c [n], defined by, 

E c [n\ = (*|f + V ee \^) - ^> s \f + V ee \^ s ) (36) 

where * yields n and minimizes the expectation value of (T + V ee ). The scaling of this 
functional is more complicated. Write 

¥,(n, . . . , v N ) = e 3 ^ 2 $(r ie l9 , . . . , v N e w ) (37) 

where $(r 1; . . . , r N ) is unknown but yields n(r). Using what we have already learned above 
about the scaling of the kinetic and Hartree functionals, 

(* 9 |f + V e e\^e) = e 2l9 ($|f + e~ w V ee \$) (38) 

If $ is then chosen to minimize the real and imaginary parts of the bi-expectation value of 

f + e~ i6 V eei 

E c [n e (r)} = (* e \f + V ee \^ e ) - {%\f + V ee \%) 
= e 2W mf + e- ie V ee \§) 

-e 2W ^ \f + e' w V ee \^o) 

= e 2i6 E c M (39) 
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where A = e , and E G: \ is the correlation energy functional for a coupling-constant Hamil- 
tonian: 

H(\) =f + XV ee + J dr v(r, A)n(r) (40) 

and v(r,\) is chosen such that the density of the ground-state is exactly n(r) independent 
of A. The complex correlation function of DFRT would then obey 

E e c [ne\ = e- 2i6 E c , x [n e \ (41) 

Note, it is not proven that the potential v(r, A), and hence the coupling-constant system, 
exists for this complex choice of A. This system is quite different from the coupling-constant 
system typically used to derive a scaling relationship for the ground-state DFT correlation 
functional. One can see that H(\) is non-Hermitian and therefore the bi-expectation value 
is used in Eq. 39. v (r, A) must be chosen such that the complex-valued function $ minimizes 
the real and imaginary parts of the bi-expectation value of T + e~' t9 V ee and yields n(r), a 
real density, from the bi-expectation value ($|n(r)|$). Assuming the existence of such a 
system shows a possible mapping between the complex correlation functional of DFRT and 
the correlation functional of ground-state DFT. Therefore, the natural choice of ground- 
state functionals evaluated at the complex density for DFRT seems reasonable provided the 
correct scaling factors are included. 

Even though the use of ground-state functionals in a DFRT implementation seems promis- 
ing, the analytic form of some approximate functionals could create issues. For example, 
consider the simple LDA appoximation to the exchange energy, 

E^[n e ] = C x J dvnfiv) (42) 

where C x = — (3/4)(2/7r) 1 / 3 . The (4/3) power in the functional is not uniquely defined for 
complex functions. This problem was considered recently by Ernzerhof et al. who concluded 
that the non-uniqueness is removed if the complex density behaves continuously as the 
imaginary part of the potential is reduced to zero 38 . One should be careful to analyze the 
form of approximate ground-state functionals before their implementation in an approximate 
DFRT. 

Along with possible uniqueness issues, any approximate functional used in an implemen- 
tation of DFRT could render a bound system unbound or an unbound system bound. In 
other words, errors in the approximate DFRT functionals could give the energy of a bound 
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state a non-zero imaginary part, hence indicating that the state has a finite lifetime. In this 
case, the imaginary part of the energy measures the error induced by an approximate com- 
plex functional. These types of errors are expected to be most visible when the ground-state 
of a system is weakly bound or for a metastable system that has a very long lifetime. 



VII. CONCLUSION 



In this work, the fundamental quantities of DFRT and practical considerations of its 
implementation have been explored. Results of model calculations show that sensitivity 
to the complex-scaling parameter can be avoided via the use of a non-interacting auxilary 
system, the Kohn-Sham system. These one-particle equations can be liberated of any 6- 
dependence by using a large, but still reasonable, basis set or fine grid. The stucture and 
significance of the orbital energies and lifetimes of this Kohn-Sham system have been studied 
for some relevant cases. Some systems, such as the LER of a system with a NEA, might have 
a mix of both bound and metastable Kohn-Sham particles. The HOMO energy is related 
to the NEA, the width of the LER, and the threshold energy of the system. In addition, 
physical interpretation has been assigned to features of complex density functions such as 
oscillations which indicate geometrical bias of decay, and a reactivity theory based on these 
functions has been motivated via model system calculations. Lastly, the complex-scaling of 
ground-state DFT functionals has been studied, a neccessary step towards general use of 
DFRT. Scaling relationships have been derived which show how the DFRT functionals can 
be written as the DFT functional evaluated at the complex density multiplied by a complex 
constant. An implementation of DFRT for real systems of interest and a time-dependent 
extension are forthcoming. 
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